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We consider a simple model, based on currently accepted models for active galactic nu- 
■ clei, for a quasi-stellar object (QSO or "quasar" ) and examine the influence that correlation- 

induced spectral changes ( "The Wolf Effect" ) may have upon the redshifts of the optical 
emission lines. 
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The Wolf effect is the name given to several closely related phenomena in radiation physics 
dealing with the modification of the power spectrum of a radiated field due to spatial fluctuations 
of the source of radiation. The first paper to demonstrate a connection between spatial source 
fluctuations (as characterized by the cross-spectral density function) appeared in 1987 |lj] . Since 
then over 100 papers have appeared on this and related topics (for a recent review, see ref. [[|). 

The discovery of coherence-induced spectral shifts grew out of work investigating the con- 
nection between optical coherence theory and the empirical laws of radiometry, in particular the 
dependence of radiative transfer on the wavelength of light. In 1986, Wolf introduced the scaling 
law |3J, which is a condition upon the spatial coherence properties of a source under which the 
radiated field has the same normalized spectrum as the source; in other words, if the source 
obeys the scaling law, the radiation pattern is independent of wavelength. From this discovery, 
it was a comparatively short step to investigate the properties of a radiation source which does 
not obey the scaling law; it was found that for such a source there would be a small shift of 
the central frequencies of spectral lines, with shifts towards lower frequencies predominating M. 
This effect was soon confirmed experimentally ||, f|, [5], ||. 

Because of the well known analogy between radiation and scattering, it is not surprising that 
analogous effects may occur in the scattering of light. Such spectral shifts were first investigated 
in ref. [0j for the case of static (i.e. time invariant) scattering media; the analysis of the case 
of time-dependent scatterers followed soon after ||. This led to an important realization that 
scattering from time-dependent fluctuating random media could mimic the Doppler effect in 
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many of its important features ||. An example of such a scattering medium which could produce 
Doppler-like shifts was found in 1990: it involved a highly statistically anisotropic [] scattering 
medium in the sense that the coherence length of the fluctuations of the scatterer was much 
longer in one direction than in the others [jl(| (see also |11 , [T2|, [13| ) . 

It was speculated in ref . that this effect might play a role in the interpretation of the spectra 
of certain extra-galactic astronomical objects such as quasars. Since then other applications, 
notably remote sensing, communications and filtering have since be proposed; However in this 
paper we will return to this first application. In the work on scattering from anisotropic media, 
it was pointed out that anisotropies of the type under consideration are compatible with current 



models of quasars |10| . Here we will investigate these points further by proposing a model for a 
quasar, closely related to a model currently favored |14]], which predicts a sizable spectral shift 
due to the Wolf effect. The rest of the paper is organized as follows: in section 2 we discuss 
briefly various types of spectral shift phenomena and their importance to cosmology. In section 
3 we give a qualitative discussion of the model quasar we are considering and a mathematical 
analysis of the spectrum of the light it radiates. Finally we assess the possible implications of 
the results derived. 



2 Spectral shifts and Cosmology 

In this section we will give a brief and non-exhaustive description of different physical phenomena 
which can give rise to shifts in the central frequency of a spectral line. The magnitude of such a 
shift is often measured in terms of the fractional shift or ^-number, defined by the formula 

Aq~-Aq CJq - uj' 

where Ao is the unshifted central wavelength of the spectral line; A is its shifted central wave- 
length; cuo is the angular frequency corresponding to Ao and u' is the angular frequency cor- 
responding to A . When 2 is a positive number the spectral line has been shifted to a lower 
frequency (i.e. a redshift) and when z is negative the shift is to a higher frequency (a blueshift). 
From eq.( |2.1| ) it is clear that z must always be in the range — 1 < z < 00. 

Possibly the most famous physical phenomenon which can produce spectral shifts is the 
Doppler effect. When a source of radiation is in motion with respect to the observer, there is a 
shift of the wavelength of spectral lines with respect to the wavelength that would be observed if 
the source were at rest. In this case the z-number is given by the the following formula Wi 



P C ° 89 l«-/3cos*, (2.2) 



where (3 — v/c,v being the speed of the source relative to the observer and c the speed of light, and 
6 is the angle between the velocity of the source and the line joining the source and the observer. 
The approximate expression is valid when /3 -C 1. For example, if the source is moving directly 
away from the observer, so that 9 = tt, the ^-number has the property z > 0, i.e. a redshift is 
observed, whereas if the source is heading directly toward the observer (so that 6 = 0) then the 
^-number is negative and a blueshift will be observed. An important property of the Doppler 



*we will use the term statistically anisotropic for such media to avoid confusion with other types of anisotropies, 
for example media which respond to different polarizations of light in different ways 
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effect is that the z-number is independent of the wavelength of the line in question. The fact that 
redshifts with ^-numbers approximately independent of wavelength are observed in the spectra 
of all galaxies, thereby suggesting that the spectral shift is due to the Doppler effect associated 
with receding motion, is the single most compelling piece of evidence in favor of the expanding 
universe. Hubble's law, the empirical linear relationship between distances and redshifts of close- 
by galaxies, confirms the predictions of cosmological models based on general relativity, and so it 
has become widely accepted that the redshift of the radiation emitted by an object is an indicator 
of its cosmological distance. However there is observational evidence that seems to contradict 
these assumptions JH]]. A well known example is the close pairing of the galaxy NGC-4319 (with 
redshift z = 0.006) and the quasar Mk-205 (z = 0.07), which, despite having drastically different 
redshifts (and therefore being a completely different distances according to the usually accepted 
interpretation of redshifts), have a distinct luminous connection between them. If we abandon 
this conventional interpretation, we must then ask what phenomena, other than the recessional 
motion of an object due to cosmological expansion, may be the cause of the observed spectral 
shifts. 

Another important type of spectral shift is the gravitational shift discovered by Einstein |TS| . 



This occurs when there is a difference between the gravitational potential in the vicinity of the 
source and the gravitational potential in the vicinity of the observer. In this case the z-number 
is given by 

z = — 2.3 

c z 

where A$ is the difference in gravitational potential. This shift, like the Doppler shift, has a 
z-number which is independent of wavelength. However, it can be shown that this shift is small 
for quasars. 

Another type of shift, which is sometimes suggested as a possible source of spectral shifts 
in quasars, is the Compton effect, due to scattering by free electrons [|19|, [20f] . The scattered 
radiation has a spectral shift with respect to the incident radiation given by the formula: 

z = 2^q sin 2 (9/2) 

m e c 2 

where m e is the mass of an electron and is the scattering angle. However, shifts due to the 
Compton effect have ^-numbers which are dependent on the wavelength of incident radiation 



(i.e. u! , the angular frequency of the the incident radiation, appears in eq. [2~4|) , and therefore the 
Compton effect cannot account for all of the large shifts observed in the spectra of extra-galactic 
objects. 

Other important optical effects which produce spectral shifts are the closely related scattering 
phenomena known as Brillouin scattering ]3T[ and Raman scattering [22, |23|. Both of these effects 



produce changes in the wavelength of light due to interactions with excitations of internal degrees 
of freedom of a scattering medium, such as rotational excitations of molecules or vibrations of a 
crystal. However these effects do not produce simple spectral shifts: the spectrum of scattered 
monochromatic radiation contains various discrete frequencies (a triplet in the case of Brillouin 
scattering, multiplets in the case of Raman scattering). 

One other hypothesized type of spectral shift should also be mentioned, namely the redshift 



due to fundamental particles having variable masses [p4| , p5| , p6[ . Although the effect has never 
been confirmed experimentally, it can be put on a sound theoretical basis, and it does seem to 
offer a possible explanation for many otherwise contradictory observations of the universe [O . 
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To these possibilities, we now add the Wolf effect, already discussed briefly in the introduction. 
In the following section we discuss a model of a quasar which displays a redshift of its spectral 
lines due to this effect. 



3 A model for active galaxies with discordant redshifts 

Quasars are star-like objects which have many unusual properties. They emit a non-thermal 
continuum spectrum of radiation and emission lines that are highly redshifted (z-numbers of up 
to 4 or 5 have been reported). Because of this high redshift, they are conventionally thought to 
be very distant objects, which therefore can reveal important information about the structure of 
the universe on very large scales. However, given their observed intensity, if this interpretation of 
their redshifts were correct, then quasars must be very powerful sources of radiation indeed. Since 
the discovery of quasars over 30 years ago, very elegant theories have evolved to explain their 
properties and account for this peculiar intensity of radiation. The possibility that the observed 
intensity might be the correct indicator of distance (rather than redshift), and that therefore 
these object might have a peculiar non-cosmological redshift is usually implicitly rejected. 

The basic features of currently fashionable models of quasars ||14j| , and other types of peculiar 
extra-galactic object (collectively known as active galactic nuclei or AGNs), are shown in figure 
1. The most important point is that these objects have axial rather than spherical symmetry. At 
the center is a very dense object (the 'central engine', sometimes characterized as a black hole, 
although no direct confirmation of the predictions of general relativity have been seen) which, 
due to gravitational attraction, sucks in matter from the dust torus, causing it to radiate with a 
broad, non-thermal power spectrum. This radiation excites gaseous clouds which will emit line 
radiation. Usually there are thought to be two types of clouds, the broad line regions (BLR) and 
the narrow line regions (NLR), which emit spectral lines of different widths. (In our diagram we 
have lumped the BLR and NLR together). The beauty of this model is that the properties of 
different types of AGN can be explained quite simply by the fact that one is observing the same 
type of object from different directions. The dust torus blocks out certain wavelengths in certain 
directions. 

Our model differs from that usually considered in two important ways. First, the line emitting 
clouds are considerably closer to the central engine (or alternatively the dust torus has a consid- 
erably larger bulge), so that the line radiation is collimated as well as the continuum radiation. 
Secondly we assume the existence of a time fluctuating scattering medium, whose properties 
will be discussed below, which causes the line radiation to be scattered into the direction of the 
observer. A schematic version of this model, suitable for mathematical analysis, is shown in fig. 2 

3.1 spectrum of scattered radiation 

Scattering by a time-dependent medium is described by the following inhomogeneous partial 
differential equation: 



In this formula, U(r,t) is the radiation field, treated in the scalar approximation, and t](r,t) is 
the time-dependent dielectric susceptibility of the scattering medium. We have used a simplified 
model for the medium, in which rj(r, t) has only a single time argument; more generally, it should 
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have two time arguments in order to account for the possibilities of resonant frequencies of 
the scatterer. 

We will analyze the scattering in the space-frequency domain. The Fourier transforms f\ of 
U(r,t) and r)(r,t) are defined by the formulas: 

1 r°° 

U(r,u) = — / U(r, t) exp(iut)dt, 

1 f°° 

r](r,u) = — / i](r, t) exp(iut)dt. (3.2) 

2ll J-oo 

Substituting these definitions into eq. (|3.1|) we obtain the following partial differential equation: 

(V 2 + k 2 ) U(r, a) = -4tt f°° fj(r, u - u')U(t, u')du', (3.3) 
where k = u/c. This can be solved using the standard Green's function: 

+ J J d 3 r'du'ri{r', u - u')U{r', J) ~ ^ , (3.4) 

where Uo(r,uj) is the incident field. This is a Fredholm integral equation of the second kind. It 
can be solved using the standard Born series. Retaining the first two terms only, we obtain the 
following expression for the scattered field: 



U s (r, u) = U (r, lu) - U (r, u) 

' rf 3 rWr/(r', u - u')U {r', u > f x ^ k ' r ~ r 'D , (3.5) 



In the far zone (i.e. in the limit |r| 3> |r'|), the scattered field is given by the formula: 
U^(ru,uj) « !^P(!M f f d 3 r'duj'fj{r',uj~uj')U {r',uj')exp{-tku-r'). 



(3.6) 

where u is the unit vector in the direction of observation (fig. 2). The power spectrum of the 
far-zone scattered field is therefore: 

S {co \ru,u) = — J [d 3 r[d 3 r' 2 j du'W^v'^u - u') 

xWoK,r' 2 X)exp[ifcu- (r[ - r' 2 )] (3.7) 

where we have used the following identities involving the cross-spectral density functions for the 
incident field, Wo(r[, r' 2 , uj) , and for the dielectric fluctuations W n {r' x ,Y 2 ,uj) and the spectrum of 
the scattered field S^°°'(ru,uj): 

(i/ s (oo) *(ru,w)i/ s (M) (ra,w')) = S ioo) {ru,ui)6{u-Lo') 

(tf Vi,^o(r>')) = Wo&^uWu-u') 

(fj*(T f 1 ,u)ij(T f 2 ,u')) = Wn&^uMu-u'), (3.8) 

the chevron brackets denoting an average over ensembles of both the field and dielectric fluctua- 
tions. 



2 Since both U(r,t) and rj(r,t) are stationary random processes, generalized harmonic analysis is required. 
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3.2 The Incident Field 

Assume that the incident field is radiated by a three dimensional primary source distribution 
Q(r,t). Further, assume that the scattering region is in the far zone of the source. The cross- 
spectral density of the field is then given by 



H/(r lSl ,r 2 s 2 ,u;) = J J ^^(f^) 



exp[zA;(r 2 — ri) 
r x r 2 

xexp[ijfe(s 1 -rj-s 2 -r£)]> (3.9) 



where the cross-spectral density function Wq{y' 1 , r 2 , u) of the source distribution is defined by 

[Q%r[,u)Q(r' 2 ,u')) =W Q (r[y 2 ,u;)d(u; - u') (3.10) 



where Q(t',u) is the generalized Fourier transform of the source distribution Q(r,t). 

It will be convenient for the later development to express this result in terms of "sum-and- 
difference" coordinates, as follows: 



x exp 



(3.11) 



where 



and 



)V(R.W) = IT (R-~b,R+-b,w) (3.12) 



r x si + r 2 s 2 

R = its = ; b = r 2 s 2 - rjSi 

R" = ^±^; r" = r£-r£. (3.13) 

If we assume that the source is quasi-homogeneous and has a normalized spectrum that is inde- 
pendent of position, i.e. 

W Q (R", r", u) = v coh s (uj)I Q (R") 5{r") (3.14) 

where Iq (R") is the intensity distribution, v co h is the coherence volume (which is formally infi- 
nite) and sq(u>) is the normalized spectrum then eq.( |3.11|) can be re- written as follows 



,.. /n . s exp(i/cs-b) 



x / d 3 ii!"J Q (R") exp ( ikh ^ " ^ ) (3.15) 



If we assume that the source is a spherically symmetric Gaussian function, viz. 



1/2 N 



/ Di 

7 Q (R'')=I exp -— 1, (3.16) 
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where a can be considered to be the effective radius of the source, then we obtain the following 
expression for the cross-spectral density function of the field in the far zone of the source: 

W(ife,b, W ) « ^^ So Mexp^s-b)exp (~°^\ (3.17) 

where Vq = (2n) 3 ^ 2 a 3 is the volume of the source. We will treat the baffle shown in fig.2, which 
represents the blocking effect of the dust torus, with a purely geometric approximation. Thus 
the cross-spectral density function for the field incident upon the scatterer is given by 

j Vcoh ¥? lQ s Q {u) exp(i£;s • b) exp , lit region 

W„(ife,b,u;) = I \ 2R J (3.18) 

[ 0, shadow region. 

3.3 The Spectral Scattering Kernel 

We will now introduce a model to describe the correlation properties of the scattering medium. 
Let us assume that we can approximate the cross spectral density of the scatterer by the following 
formula: 

W„(riy 2 ,o;)«p(i2)§(b,s,a;), (3.19) 

where 

2 

b = r'a-ri. (3.20) 

This model is somewhat analogous to the well-know quasi-homogeneous model, with the added 
wrinkle that the coherence properties can be dependent on orientation (as characterized by the 
unit vector s). What I have in mind specifically is that the coherence length in the radial direction 
(i.e. along s) is different from the coherence length in the transverse direction (i.e. perpendicular 
to s). 

Therefore the spectrum is 

S (oo) (ru,a;) = —j J d 3 Rd 3 b J dw'p(R)g(b, s, cj)W (R, b, u') exp(zA;u ■ b). (3.21) 



On substituting from eq.( |3.18j ) into eq.( |3.21| ) we obtain the following expression for the spectrum 
of the scattered radiation 



S (oo \ru, uo) = —du f K(u, u, u/)s Q (u/) (3.22) 

where the spectral scattering kernel k(u, ujuj') is given by the formula 

r r ( k' 2 a 2 b 2 \ 

/ d 2 tt s / d 3 b g(h, s, u - u') exp exp [i (k's - fcu) ■ b] (3.23) 



Kq 

where k! = uo' / c and 



«o = VcohVoIo / dRp(R). (3.24) 
Jr 

In these formulas, Rq represents the distance from the source to the baffle and Qq is the solid 
angle of the lit region (see fig.2). 
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3.4 The Scattering Medium 

We shall assume that the scattering medium is anisotropic, in the sense that the correlation 
length in the radial direction (i.e. along the vector s) is different from the correlation in the 
transverse direction (i.e. perpendicular to s). For simplicity, we will assume a Gaussian model, 
i.e. 

' h 2 bl 



g(b, r) = exp 



en 



T 2 C 2 



(3.25) 



where = (s.b)s and bii = b — b^. Substituting this formula into eq. ( |3.23| ), we obtain the 
following expression for the scattering kernel 



■ exp 



x 



/ d 2 Q s f d 3 b exp 

J Sin J 



1 ( >\ 2 




exp (if • b) 



(3.26) 
(3.27) 



where 
and 



f = k's — ku 



a 



12 



1 a 2 k' 2 
— + 



R 2 



en 



(3.28) 



(3.29) 



The approximation made in the last formula is valid provided that the traverse coherence area of 
the incident light is much larger than that of the scattering medium. Given that the propagation 
distance R is of the order on tens or even hundreds of parsecs, this is not an unreasonable 
assumption. Performing the integration over b we obtain the following formula for k: 



k(u, u), uj') 



2ikj t <j 2 l <t\\ 
k exp 



lot 



1 r 



x I d 2 Q s exp ( — - a]k u]_ + a 2 (k' - ku\\ 

(27r) 2 oy<7i<7|i 
k 1 exp 

J cos 9q V 2 

+ a\{k' - kxcosO') 2 



1 <*T < >v 



aik z (l-x z cos 2 9') 



(3.30) 



where 6' is the scattering angle, as shown in fig. 2. The integral in eq.( |3.30| ) can be evaluated 
approximately using the following formula: 



cos 8q 



f(x)dx 



(l-COS0o)/ 

— f ( I — — 

2^r V 4^ 



1 + cos 6n 



(3.31) 
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Then we obtain the following 



k(u, u, lo ) ~ No- 



where 



a = — 



P = To 



a 



exp - 



olJ 2 - 2(3uu' + au £ 



(3.32) 



2 , 2,2/2 2 



2 , 2 



^ 2 | _2 



1 ^<M /J' 
1 ) COS# 

4:71, 



(3.33) 



3.5 The z-number of the scattered radiation 



The scattering kernel given in eq.( p.32[ ) is equivalent to one we have already investigated in a 
previous paper [see ref.]12|], eq.(12)]. In that paper, we showed that such a spectral scattering 
kernel causes Doppler-like spectral shifts with ^-numbers given by the formula: 



a 

W\ 



- 1. 



Substituting the expressions for a and (3 from eq.( |3.33| ), we obtain the following result: 



1-g 



(3.34) 



(3.35) 



Let us assume a 'white noise limit' for the scattering medium: i.e. its time fluctuations 
are effectively uncorrelated: a T — > 0. Then the formula for the spectral shift of the scattered 
radiation is: 



z(e,q) = - + (l-e)q 



Q 



where 




1 ) cos . 

Ati, 



(3.36) 

(3.37) 
(3.38) 



Equation Q3.361 ), which is displayed graphically in Fig. 3, is our final result. It demonstrates that, 
for the specific model of active galactic nucleus we are considering, the spectral shift observed will 
be dependent on the viewing angle 9', the solid angle of the radiation cone Q and the statistical 
anisotropy of the scatterer, as represented by the parameter e. When e = 1, the scattering 
medium is statistically isotropic: the coherence lengths are the same in all directions. In the 
limit e — >• the longitudinal coherence length is much greater than the transverse coherence 
length, i.e. one has a cigar shaped coherence volume. In the limit e — > oo, one has the opposite 
case: the longitudinal coherence length is much smaller than the transverse coherence length and 
the coherence volume is like a pancake. Note that z(e, q) > provided that e > 0.5, so that one 
will have a preference redshifts over a large area of the (e, q) plane. 
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4 Discussion and Conclusions 



One can easily spot potential problems with the theory presented here. For example, we have been 
deliberately vague about the underlying physical nature of the scattering medium, other than 
specifying its anisotropic coherence properties. The scatterer, which is assumed to have a 'white 
noise' power spectrum (implying that its fluctuations are very energetic), is situated further away 
from the central engine than the line emitting clouds (which cannot be too hot, otherwise they 
would have completely ionized, making line radiation impossible). Although our results apply 
for a broad spectral range, we have not considered the shifts of absorption lines or of the 21cm 
radio line. We have ignored the unscattered radiation and the efficiency of the scattering process 
(although it is possible a stimulated version of the spontaneous scattering process considered 
here is applicable). Also the issue of spectral linewidths has not been addressed. However, our 
results do demonstrate the plausibility of the existence of extra galactic objects with discordant 
redshifts due to well established physical phenomena. 

It should be emphasized that our theory is based entirely on established principles of optics: 
no hypothetical new physics was introduced. Furthermore certain aspects of the underlying 
phenomenon, the Wolf effect, have been extensively tested in the laboratory. The formula (|3.36|) 
gives the relationship between the z-number of the observed radiation and certain variables 
dependent on the object, namely the angle of vision, the degree of collimation due to the dust 
torus, and the degree of anisotropy of the scatterer. As it may be possible to measure some of 
these quantities for certain objects, it therefore may be possible to test these predictions. 
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Figure Captions 

Fig.l: A schematic illustration of the model AGN under consideration. 

Fig. 2: A simplified version of Fig.l, for mathematical analysis, showing the symbols used. 

Fig. 3: The relationship between z number and the collimation factor q for different anisotropy 
factors e, illustrating the result eq. (|3.36| ). 
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